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ABSTRACT 

We report on two regularly rotating galaxies at redshift z ^ 2, using high-resolution spectra of the 
bright [Cn] 158/4111 emission line from the HIFI instrument on the Herschel Space Observatory. Both 
SDSS090122.37+1 8 1432.3 (“S0901”) and SDSSJ120602.09+5 14229.5 (“the Clone”) are strongly lensed and 
show the double-horned line profile that is typical of rotating gas disks. Using a parametric disk model to fit 
the emission line profiles, we find that S0901 has a rotation speed of usin(f) 120 ± 7kms _1 and a gas ve- 
locity dispersion of <r g < 23kms _1 (ler). The best-fitting model for the Clone is a rotationally supported disk 
having usin(i) ~ 79 ± I 1 km s~ 1 and 4kms 1 (la). Flowever, the Clone is also consistent with a family 
of dispersion-dominated models having a g = 92 ± 20 km s ' 1 . Our results showcase the potential of the [C n] line 
as a kinematic probe of high-redshift galaxy dynamics: [Cn] is bright, accessible to heterodyne receivers with 
exquisite velocity resolution, and traces dense star-forming interstellar gas. Future [C ii] line observations with 
ALMA would offer the further advantage of spatial resolution, allowing a clearer separation between rotation and 
velocity dispersion. 
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1. INTRODUCTION 

Typical star-forming galaxies in the nearby universe are 
supported by systematic rotation with random velocities playing 
a relatively minor role. This indicates a degree of maturity — the 
galaxies are old and established, and their reservoirs of star- 
forming gas are dominated by material that was accreted at 
least one orbital time ago. However, at the epoch when cosmic 
star formation reached its peak, it is likely that galaxies were 
still accreting gas rapidly. The relative roles of rotation and 
velocity dispersion in providing support to such galaxies could 
thus be significantly different. Certainly, the morphologies of 
star-forming galaxies at z ~ 2 suggest some differences in their 
typical properties: regular spirals constituted a much smaller 
fraction of galaxies at z ~ 2 than they do today. 

Dynamical observations of star-forming galaxies at z ~ 2 are 
challenging. One successful approach has been to study strong 
rest-frame optical emission lines, which are redshifted to near- 
IR wavelengths (e.g., Forster Schreiber et al. 2009; Lehnert et al. 
2013; Rhoads et al. 2014). Spatial resolution is generally limited 
by seeing to >0 5 , which yields a few resolution elements across 
the disk of the galaxy. For a minority of cases it has been possible 
to achieve significantly higher spatial resolution, either by using 

* Herschel is an ESA space observatory with science instruments provided by 
European-led Principal Investigator consortia and with important participation 
from NASA. 


adaptive optics (e.g., Law et al. 2009), or by observing strongly 
lensed galaxies (Jones et al. 2010; Frye et al. 2012; Wuyts et al. 
2012; Jones et al. 2013). Spectral resolution is also usually 
modest ( R < 6000, corresponding to Av > 50kms _1 ) since 
the faint flux levels of z ~ 2 emission lines usually preclude the 
lower sensitivity of high-resolution spectrographs. 

Here we adopt a complementary approach to studying the 
dynamical support of z ~ 2 star-forming galaxies. We use [C n] 
158 /im line observations from the Herschel Extreme Lensing 
Line Observations (HELLO) program (Malhotra et al. 2014), 
obtained using the HIFI instrument (de Graauw et al. 2010) 
on the Herschel Space Observatory (Pilbratt et al. 2010). The 
resulting line profiles are at exquisite spectral resolution, though 
entirely unresolved spatially. In the present paper we analyze 
the galaxies SDSS090122.37+181432.3 (Diehl et al. 2009, 
hereafter S0901) and SDSSJ120602.09+5 14229.5, also called 
“the Clone” (Lin et al. 2009). Both galaxies show the double- 
horned rotation profile that is characteristic of rotationally 
supported disks. The shape of this profile contains considerable 
information on the disks’ dynamical properties. 

C + is the main coolant of neutral gas in galaxies. Typically, 
0.3% of far-infrared luminosity emerges in the [Cii] 158 /im 
line. Thus, it can be an excellent tracer of interstellar medium 
(ISM) kinematics. The question then is where the [C n] emission 
comes from. [Cii] can trace neutral gas, since carbon has a 
lower ionization potential than hydrogen. The brightest [Cii] 
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emission in a galaxy comes from photon-dominated regions 
(PDRs). These are surface layers of dense molecular clouds, 
where UV light from nearby star formation is dissociating 
molecular gas. [Cn] probes primarily a density range from 
a few cm -3 , which is required to excite the transition, up 
to the critical density of ~10 35 cm' 3 . The integrated [Cn] 
emission from whole galaxies includes potentially important 
contributions from diffuse ionized gas (Bennett et al. 1994) and 
also from diffuse Hi, although the fractions are debated (e.g., 
Madden et al. 1993; Heiles 1994; Contursi et al. 2002). The 
radial distribution of [C ii] in nearby galaxies generally follows 
CO (e.g., Kramer et al. 2013; Pineda et al. 2013, forM33 and the 
Milky Way, respectively), although [C n] emission has also been 
observed in outflows (Contursi et al. 2013) and in turbulently 
heated gas in nearby radio galaxies (P. Guillard et al. 2014, 
in preparation). The galaxies studied in this paper have typical 
[Cii]/FIR ratios (Malhotra et al. 2014), and do not show the 
[Cii] deficiency seen in some luminous nearby systems (e.g., 
Malhotra et al. 1997, 2001; Carilli & Walter 2013). 

We describe our sample and observations in Section 2. We 
discuss our Markov chain Monte Carlo (MCMC) spectral line 
fitting methods and results in Section 3, and compare them 
to Ha observations of the same objects in Section 4. We 
close with a discussion of implications and future directions 
in Section 5. Throughout the paper, we adopt a A-CDM 
“concordance cosmology” with CIm — 0.27, il.\ — 0.73, and 
H 0 = 71 kms -1 Mpc -1 . 

2. THE SAMPLE AND THE OBSERVATIONS 
2.1. Sample Selection and Properties 

The full HELLO sample consists of 15 strongly lensed galax- 
ies, at redshifts 1.0 < z < 3.0. The sample was selected to 
span a range of both star formation rate (SFR) and redshift. All 
HELLO targets were observed using the HIFI Wide Band Spec- 
trometer (WBS) backend. We describe the HIFI observations of 
the full sample in Malhotra et al. (2014). Here, we focus our 
attention on the two objects (S0901 and the Clone) that show 
distinct double peaks in their line profiles. We summarize the 
key observational information for these two objects below. Both 
were identified as candidate lens systems through a search for 
blue objects near luminous red galaxies (LRGs) in Sloan Digital 
Sky Survey (SDSS) imaging, followed by visual inspection for 
arc-like morphology (Diehl et al. 2009; Lin et al. 2009). 

S0901 was confirmed as a strongly lensed galaxy at redshift 
Z = 2.2558 by Diehl et al. (2009). The lens is a galaxy group, 
dominated by a LRG at z = 0.3459 (Diehl et al. 2009), and 
containing another half dozen smaller galaxies. The lensed 
galaxy appears as a 1 7" long arc, oriented primarily north-south, 
and located east of the lensing galaxy. Diehl et al. (2009) identify 
two main components of the arc (labeled “b” and “d” in their 
published image), and measure its Einstein radius as 1". 1 ± 1" 1. 
Hubble Space Telescope ( HST) imaging (from program 1 1602, 
PI: S. Allam) shows considerable substructure (as is usual 
for star-forming galaxies observed in the rest-frame UV), as 
well as a possible counterimage 4" west of the LRG, and a 
foreground galaxy superposed on arc component “b.” Hainline 
et al. (2009) published the lensed galaxy’s rest-frame optical 
spectrum, which shows prominent emission lines of Ha and 
[N ii], as well as weaker lines of [O in], H/8, [O n], and [Ne in]. 
Hainline et al. (2009) note that [N n]/Ha > 0.65, implying the 
presence of an active galactic nucleus (AGN). On the other hand, 
the object’s mid-infrared spectrum shows strong polycyclic 


aromatic hydrocarbon (PAH) emission features, suggesting that 
the dust heating is dominated by star formation rather than 
by the AGN (Fadely et al. 2010). The gravitational lensing 
amplification of S090 1 does not appear to be formally published. 
Available estimates range from amplification /x ~ 8 (Fadely 
et al. (2010), quoting from work in progress by A. West et al.), 
to /x ~ 6 (E. Buckley-Geer 2013, private communication), 
while inspection of the HST imaging suggests that larger 
amplifications are plausible. Where needed, we use a fiducial 
value /x = 6, but we retain /x in our equations to show explicitly 
how the results will scale for other /x. Apart from the inferred 
SFR of S0901, our primary conclusions in the present paper are 
essentially independent of the lensing amplification. 

The Clone (SDSSJ120602.09+514229.5) was first identified 
as a strongly lensed galaxy at redshift z = 2.00 by Lin et al. 
(2009). The lens is again a group, here with three identified 
members, dominated by an LRG at z = 0.422 (Lin et al. 
2009). The arc has three primary components which together 
span about 120°, plus a small counterimage. Lin et al. (2009) 
labeled the components of the main arc as Al, A2, A3, and the 
counterimage as A4. They modeled the lensing geometry based 
on ground-based imaging, finding an Einstein radius 3" 82, and 
an amplification factor /x = 27 ± 1. The lensed galaxy’s rest- 
frame optical spectrum (Hainline et al. 2009) is that of a star- 
forming galaxy, with no evidence for AGN activity. Detected 
lines include Ha, H/l, [On], [Om], and [Sii]. The Clone’s 
mid-IR spectrum is again that of a normal star-forming galaxy, 
with prominent PAH emission features (Fadely et al. 2010). 
More recent lens modeling, using higher-resolution HST images 
(Jones etal. 2010), gives a refined amplification /x = 28.1 ±1.4. 
(This is the total amplification for all lensed images combined, 
as appropriate for our spatially unresolved HIFI observation.) 

Hainline et al. (2009) also estimated the SFR for the Clone. 
They scaled from the Ha flux, assumed a Chabrier (2003) 
initial mass function, and corrected for gravitational lensing 
amplification (/x, see below). They obtained SFR = 32 M 0 yr 1 
before dust correction and SFR = 64-75 Mq yr - 1 after dust 
correction (depending on the adopted reddening indicator). They 
did not estimate a SFR for S0901, given their concerns about 
the possible role of an AGN. We can obtain a crude upper limit 
on the SFR of S0901 by neglecting AGN contributions to Ha 
and applying the same scalings used for the Clone. This yields 
~135 Mq yr -1 x (6//x), prior to any dust correction. 

2.2. Herschel Observations 

We observed the HELLO sample using the HIFI WBS with 
1.10 MHz resolution. We used dual beam switching with fast 
chopping to achieve the best possible baselines in our spectra. 
For S0901, we used mixer band lb and a total observing time 
of 3000 s, of which 885 s was on-source integration time. For 
the Clone, we used mixer band lb with a total observing time of 
6200 s, of which 1 977 s was on-source integration time. Our data 
reduction — which we discuss in greater detail in Malhotra et al. 
(2014) — followed mostly standard HIFI procedures. The last 
step of this analysis was to fit and subtract a baseline from each 
spectrum. This means that the spectra used in the present paper 
are “line-only” spectra, insensitive to the presence of continuum 
emission. 

The HIFI WBS offers an intrinsic velocity resolution of 
<0.6 kms -1 at the frequency of the redshifted [Cii] line in 
our targets. However, each 1.10 MHz resolution element has 
a signal-to-noise ratio (S/N) < 1 in our data. Our ability to 
distinguish fine detail in the spectral line profiles is therefore 
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Figure 1. HIFI spectrum of the [Cn] 158 (tm line from 
SDSS090122. 37+181432. 3, binned to 25 km s -1 resolution (black his- 
togram). Plotted error bars are derived empirically from the rms of the 
spectrum, measured at line-free frequencies. The overplotted model is a 
centrally concentrated disk model, with t>d,i /2 = 121 kms -1 , a g = 29 kms -1 , 
and a negligible halo component. The resulting ratio of circular speed to 
velocity dispersion is 4.2, and at the 95% confidence level, the ratio is >2.7. 

(A color version of this figure is available in the online journal.) 

effectively limited by sensitivity, since we can detect features 
only when we combine enough frequency channels to get S /N 
of a few in the feature. 

We plot the [Cii] spectrum of S0901 in Figure 1 and of the 
Clone in Figure 2, both at 25 km s -1 binning. S0901 has a flux of 
1.3 K kms -1 (1.2 x 10~ 17 W m -2 ), and yields a S/N > 4 per bin 
for 1 1 bins above half-maximum intensity. The Clone is fainter, 
with a flux of 0.3Kkms -1 (3.2 x 10~ 17 Wm~ 2 ), and yields 
S/N ~ 3 per bin for 6 bins above half-maximum intensity. 
We discuss these flux measurements further in Malhotra et al. 
(2014), in the context of the full HELLO sample. In addition 
to the spectra, we plot model line profiles from our fitting 
procedures in Figures 1 and 2. 

3. EMISSION LINE PROFILE FITTING 

The observed [Cii] emission line profiles depend on the 
rotation curves of the lensed galaxies. They also depend on 
the spatial distribution of the emitting gas, which determines 
how the different parts of the rotation curve are weighted in 
forming the line profile. To model the line profiles, we adopt 
commonly used fitting functions for the distributions of both 
emitting material and of mass (which in turn determines the 
rotation curve). 

For the radial distribution of [C ii] surface brightness, we have 
assumed an exponential disk model (E[c H] = Zo exp(— r/rpn])), 
which is able to produce models that fit the observed lines well 
(see below). Of course, this does not mean that the gas distri- 
bution must follow an exponential distribution. Different tracers 
of galaxies’ interstellar gas often follow different distributions. 
Neutral hydrogen (Hi) tends to have a spatially extended dis- 
tribution. Molecular gas can be centrally concentrated (as in 
many starbursts and ultraluminous infrared galaxies), or can 
show ring features, as in both our own Galaxy and M31. Ha 
comes primarily from dense, ionized gas, and often shows com- 


Figure 2. Like Figure 1, but for SDSSJ120602.09+5 14229.5 (the Clone). Here, 
the signal-to-noise ratio is lower, and a range of models provide acceptable 
fits to the data. The modeling procedures are described in Section 3. We plot 
four representative cases. The dashed red line shows the model with the highest 
likelihood among our MCMC simulation results, a "pure” double-horned line 
profile where most of the emitting gas is on the flat part of the rotation curve, 
and the velocity dispersion is very small. If we consider only models with 
cr g > 7 kms -1 (motivated by typical motions within thin disks), the best fit 
(shown as the solid red line) remains double-horned, but with very extended 
wings, here driven by a rotation curve that continues rising beyond the spatial 
extent of the emitting gas. Finally, we plot two dispersion-dominated models. 
The first (red dotted line) is effectively a pure Gaussian, with a g = 98 km s 
and negligible rotation. The second (blue solid line) matches the results of Jones 
et al. (2013), from spatially resolved near-IR spectra, having ly, 1/2 = 70 km s -1 , 
Xds =1.1, and a g = 94 km s 1 . The ratio of circular speed to velocity dispersion 
spans a wide range for these and other acceptable models of the Clone’s [C 11 ] 
emission (see Figure 3). 

(A color version of this figure is available in the online journal.) 


plex and patchy distributions that trace the distribution of young 
star-forming regions. 

A double-horned line profile basically requires that most of 
the emitting material be in significant motion with respect to the 
center of mass. Given that galactic rotation curves rise steeply 
in the central regions and then flatten out, there are two “usual” 
ways to achieve this. First, the gas can be in an extended disk, 
with most of the mass on the flat part of the rotation curve. 
This is the common case for H 1 (Springob et al. 2005). Second, 
the gas can be in a more compact ring-like configuration, with 
little material to emit in the regions where rotation is slow. 
This has been observed, e.g., for CO line emission in NGC 759 
(Wiklind et al. 1997). Since [Cii] 158 /um emission comes in 
some combination from the diffuse warm neutral medium, PDRs 
at molecular cloud surfaces, and ionized regions, we do not 
expect a simple analogy with other main ISM tracers. We begin 
by considering the extended disk possibility in detail, and then 
examine how our conclusions would be affected for the case of 
a [C 11] ring. 

We model the rotation curve using a disk + halo model for 
the gravitating mass. The halo has a mass volume density of 
p oc (1 + r/rh)~ 2 and an asymptotic rotation speed of 17,. The 
disk component has a mass surface density of Z oc exp (— r/rj) 
and a peak rotation speed of Vd . We can relate v,/ to other 
disk parameters as Vd ^ 0.5463V G M^/ rd, where Mdisk is 
the total mass of the exponential disk. Vd corresponds to the 
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circular speed for the disk mass interior to r = 1.7933rd. 12 We 
model the line-of-sight gas velocity dispersion (due to thermal 
and turbulent motions) as a Gaussian with rms width a g , which 
is convolved with the line profile produced by ordered rotation. 

Because our targets are spatially unresolved in HIFI data, our 
predicted line profile is insensitive to a global rescaling of all 
radii by a common, constant factor. We therefore fix the scale 
length r[Cnj of the surface brightness distribution throughout 
the modeling, and rewrite the scale lengths of the disk and halo 
in dimensionless form as Xd = rd/r[Cn] and x/, = r^/r^ cnj. 
Similarly, the line profile depends on the line-of-sight velocity 
v sin(/) (where i is the inclination angle of the disk) but cannot 
distinguish changes in i from changes in v. 

The line shape is then determined by a total of five free pa- 
rameters: Three velocities ( Vd , 17, , and a g ), and two scale length 
ratios (xj and a/,). We need two further parameters, the redshift 
Z and total line flux /[Cn]- to completely specify the line. 

Before fitting the disk+halo rotation curve models, we recast 
the three velocity parameters in a new form that reduces corre- 
lated uncertainties in the parameter space. Since the observed 
line profile is determined by gas kinematics where the light is 
emitted, we define the rotational speed contributions due to the 
halo and disk, Vh,i/2 and 1771/2 measured at the half-light radius 
of the [C 11] disk (which is r[c n ],t/2 = 1.66r[c U ])- We then define 
the total ordered rotation speed i) c j/2 at the half-light radius by 


We further define 

«i/2 = v l 1/2 + °g = w*. 1/2 + Vd, t/2 + • ( 2 ) 

The parameter v \/2 now contains essentially all the information 
required to determine the [C 11] line width. The relative impor- 
tance of random and ordered motions is captured by the velocity 
ratio 

ycrc ®g/Vc, 1/2 (3) 

while the relative importance of disk and halo terms in the 
potential is captured by 

ydh — Vd,\/ 2 /Vh,l/ 2 ■ (4) 

The transformation is essentially from Cartesian coordinates 
(u/,,1/2, 1771/2, er) to spherical polar coordinates (121/2, y ac , ydh ) 
for the three-space of velocity parameters, with the y parameters 
corresponding to (co)tangents of two angles. 

The final set of parameters for fitting then becomes 
{vi/ 2 ,y<rc,ydh,x d ,x h ,z, /[cn]}- By recasting the three veloc- 
ity parameters, we eliminate major degeneracies between Vd, 17, , 
and (jg, all of which have the same lowest-order effect of increas- 
ing the model line width. Now, iq/2 determines line width, while 
the remaining parameters control the shape of the line — whether 
it is singly or doubly peaked, the splitting of the peaks and depth 
of the gap between them, and the shape of the high-velocity line 
wings. 

We fitted the rotation curve models to the data using a 
MCMC approach (more specifically, the Metropolis-Hastings 
algorithm). In brief, the method consists of a chain of random 
steps in our parameter space. In each step, the algorithm first 
identifies a trial solution that is offset from the current solution 
by some randomly determined amount. The likelihood of the 

12 Rephrased, the rotation curve for a pure exponential disk peaks near 1.8 
scale lengths. 


trial solution is compared to that of the current solution. The trial 
solution is adopted as the new solution whenever its likelihood 
is higher. If its likelihood is lower, it may still be adopted as the 
new solution, with a probability equal to the ratio of likelihoods. 
After many steps, the distribution of accepted points provides 
an unbiased estimator for the probability distribution of the 
model parameters. Since consecutive steps of the chain are 
highly correlated, we record the state of the chain every 100 
steps, so that the autocorrelation among the recorded points in 
parameter space is weak. 

Over 99% of the MCMC solutions for S0901 under the 
disk+halo rotation curve model lie in a single region of pa- 
rameter space, where rotation dominates over a small or negli- 
gible amount of random motion. A representative fit is shown in 
Figure 1, and its parameters are given in Table 1. The distri- 
bution of MCMC samples in the (iq/2, er g ) plane is shown in 
Figure 3. Only 5% of the samples have o g > 43.5 kms -1 ; 1% 
have o g > 60 kms -1 ; and only 0.5% have o g > V\/ 2 - More- 
over, the rotational motions are usually dominated by the disk 
component (with (log (ydh)) = 0.92, and with ydh > 1 in 85% 
of samples). 

This conclusion is robust, being readily explained by exam- 
ination of the observed line profile. The line is double-horned 
with a central dip, and with rotation curve sides that are steep 
but not vertical. Figure 4 shows examples of the line profiles 
produced by three simple kinematic models, each of which is 
dominated by a single component. Models dominated by ve- 
locity dispersion (o g ) are expected to be single-peaked. Models 
dominated by a halo (with small core radius and without a sig- 
nificant disk contribution) would have vertical sides, except for 
the influence of the velocity dispersion. Under our adopted ro- 
tation curve model, the best way to get a double-peaked profile 
with slanting sides is to have a rotation curve that rises to a 
peak (driven by the disk component) and then declines toward 
somewhat smaller rotation speeds (perhaps dominated by the 
halo component). 

Our [C 11] line profile fitting in the Clone yields a range of 
viable models. The maximum likelihood fits are double -peaked 
models with very small velocity dispersions, much like the 
fits for S0901. However, these models occur in a rather small 
volume of the sampled parameter space. The MCMC fitting 
for the Clone also yields single-peaked, dispersion-dominated 
fits (with y ac > 1) in about 45% of all samples. These have 
appreciably lower likelihoods, but occupy a large volume in 
parameter space. For the Clone, we have plotted a few different 
models, all statistically acceptable fits to the data. These include 
the maximum likelihood model identified, which is dominated 
by a disk with 17/1/2 = 78 kms -1 and Xd s = 1.52, and has 
Og 1 km s -1 . We plot also a second double-peaked profile and 
two dispersion-dominated single peaked profiles. Parameters for 
all are summarized in Table 1. The last of these plotted models 
was selected to be consistent with the parameters that Jones 
et al. (2010) derived based on Ha observations (see Section 4). 

Our fitting yields other parameters in addition to v \/2 and o g . 
The most robust and model-independent are the redshifts (which 
we report in the local standard of rest [LSR] frame). For S0901, 
we find zjcni = 2.25860 ± 0.00007. This matches precisely the 
rest-frame optical redshift ( z = 2.2586± ~ 0.0001; Hainline 
et al. 2009) derived from Ha, [Nil], and [Om] emission 
lines, while it is slightly higher than the rest-UV redshift of 
Z = 2.2558 ± 0.0003 (Diehl et al. 2009). This gives a net blue- 
shift of Av ~ 260 ± 30 km s -1 for the UV spectral features, 
relative to [Cn]. The Diehl et al. (2009) redshift is from a 
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Figure 3. Probability distribution for the kinematic support parameters u ro t andcr^ for SDSS090122.37+1 8 1432.3 (left) and for the Clone (SDSSJ 120602.09+5 14229.5, 
right), based on our disk+halo rotation curve model (Section 3). The plotted points show nearly 8000 samples from our Markov chain Monte Carlo parameter fitting 
code. Plotted contours enclose 68%, 95%, and 98% confidence regions. For S0901, rotationally supported solutions are strongly favored, with only 0.5% of the samples 
having o g > u ro t. For the Clone, which has a lower signal-to-noise ratio, the corresponding fraction is ~44%. 

(A color version of this figure is available in the online journal.) 


Table 1 

Parameters of Selected Model Fits Plotted in the Figures 


Object 

Vd, 1/2 

I'll, 1/2 


Xd 

Xh 

Flux 

A log(L) a 

Line 


(kms -1 ) 

(kms -1 ) 

(kms *) 



(Kkms -1 ) 


Style 

SDSS 0901 

121 

2 

29 

0.11 

9.3 

1.359 

0 

Solid red 

The Clone 

78 

2 

1 

1.52 

0.7 

0.2657 

0 

Dashed red 

The Clone 

20 

100 

8 

0.07 

6.2 

0.365 

-1.4 

Solid red 

The Clone 

0 

2 

98 

0.5 

0.26 

0.371 

-3.3 

Dotted red 

The Clone 

70 

2 

94 

1.1 

0.05 

0.386 

-3.44 

Solid blue 

(None) 

75 

1 

7.5 

0.22 

0.1 



Black 

(None) 

1 

75 

7.5 

1 

0.1 



Blue 

(None) 

10 

10 

75 

1 

1 



Red 


Notes. We include the best-fitting model for S0901 (top line, plotted in Figure 1); four representative model fits 
for the Clone (lines 2-5, plotted in Figure 2); and three sample “toy models” chosen to illustrate features of 
the rotation curve model (lines 6-8, plotted in Figure 4). In each case, parameters associated with the dominant 
component(s) determining the line profile shape are given in bold, and small changes to the remaining parameters 
would have small or negligible effects on the model line profile. 

a A log(L) is the log of likelihood, relative to the maximum likelihood model identified for the same observed line. 


cross-correlation of the S0901 spectrum with a template from 
Shapley et al. (2003), making it a weighted average of the Lya 
line redshift and the redshift from UV absorption lines (which 
are, individually, mostly too faint to measure well in S0901). 
To the extent that this represents a Lya redshift, it is an unusual 
result — Lya is generally redshifted with respect to other velocity 
tracers (e.g., Shapley et al. 2003; McLinden et al. 2011), and 
a blueshift of this magnitude for Lya would be worth further 
investigation. 

For the Clone, we find z = 2.0030 ± 0.00009. This can 
be compared to the measurements z = 2.0001 ± 0.0006 and 
z = 2.0010 ± 0.0009, both from rest-UV absorption lines 
reported by Lin et al. (2009) but based on different spectra 
taken at two different telescopes. Comparing our [C ii] redshift 
to the weighted average, zuv = 2.0004 ± 0.0005, we infer that 
the gas responsible for the UV absorption lines is blueshifted 
by 260 ± 60 km s -1 with respect to the [C ii] emitting gas. This 
is typical of the absorption line blueshifts that are seen in many 


high-redshift galaxies, and that are thought to be caused by 
absorption by gas in a galactic wind (e.g., Pettini et al. 2001; 
Frye et al. 2002; Shapley et al. 2003). 

Finally, we consider how our conclusions would be affected 
by changing some of the assumptions behind our rotating disk 
model. First, we explore the possible effects of differential 
gravitational magnification by applying separate amplification 
factors to the approaching and receding portions of the galaxy. 
Taking the overall mean amplification to be jl and the ratio 
of red-side to blue-side amplification to be //„, we have 
amplifications of 2/x /( 1 +/x fl ) for the blue side and 2/r fl jl/( 1 +/x a ) 
for the red side. We apply this differential magnification to 
the predicted line profile before accounting for the velocity 
dispersion, thus avoiding a sharp step function at the systemic 
velocity. 

For S0901, we find that the distribution of this new parameter 
is sharply peaked near unity: 98.64% of simulated samples are 
in a primary peak centered at /i a = 0.92, with dispersion of 
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Figure 4. Correspondence between rotation curve and line profile for three “toy models,” each of which is dominated by a single kinematic component. In both panels, 
blue corresponds to a halo-dominated rotation curve (with gas velocity dispersion a g =7.5 km s -1 ), black to a disk-dominated rotation curve ( a g = 7.5 kms -1 ), and 
red to a dispersion-dominated galaxy (o g = 75 km s -1 ). Left: the rotation curves (with shading for the region v c (r) ± a g ). The two green curves show dL[Qn\/dr and 
its integral, L[cn](r) (i.e., the enclosed [Cn] luminosity), both normalized to their respective peaks. They show how the different portions of the rotation curve are 
weighted in the line profile. Right: the resulting line profiles. 

(A color version of this figure is available in the online journal.) 


±0.17. The 1.3% of simulations with more extreme ratios are 
all velocity-dispersion dominated models, where the asymmetry 
introduced in the magnification is hidden by the random motions 
of the gas. The only correlation between /. i a and other model 
parameters is a trivial correlation with model redshift: an offset 
in redshift that would move the modeled flux centroid away 
from the observed flux centroid can be compensated by raising 
the peak of the line profile that falls nearest the observed flux 
centroid, and lowering the amplitude of the other peak. The 
addition of does not alter the distributions of other model 
parameters in any qualitatively important way. 

For the Clone, the /x„ parameter is not strongly constrained. 
Models having a g < i> c ,i /2 generally have 0.4 < < 2.5, but 

dispersion-dominated models can yield acceptable fits with any 
value of fi a . With the addition of //„ , about 3 /4 of the fits end 
up with (7 g > v C 'i/ 2 - 

We also explored the possibility of a ring-like gas distribution. 
We take the limiting case by assuming all the emitting gas 
rotates with a single circular speed v c . The intrinsic line profile 
is then proportional to 1 /^j 1 — (v/v c ) 2 , potentially modified 
by the differential amplification parameter [i a described above, 
and convolved with a Gaussian of width a g to account for 
velocity dispersion within the ring. In total, this is a five- 
parameter model: {v c , a g , z, f[c h] , /x fl }. Qualitatively, our key 
results are unchanged under this model: we again find that 
S0901 is rotationally supported, and that the Clone admits a 
range of acceptable models. In closer detail, we do see some 
differences. The full disk models described above can smooth 
the line profile of a single ring both through velocity dispersion 
and through the combination of multiple rings with different 
circular speeds. The single-ring model can only soften the sharp 
edges of the double-horned profile through the effects of velocity 
dispersion. Thus, values of 30 km s 1 are common in the 
full disk model for S0901 but rare among the single-ring fits. 
Similarly, the distribution of y ac = a g /v c tends to slightly higher 
values under the single-ring model. Our exploration of ringlike 
gas distributions leads to two main conclusions. First, the 


Atacama Large Millimeter Array (ALMA) observations will be 
tremendously valuable in determining the spatial distribution of 
the emitting gas in cases like these. Second, our conclusion that 
S0901 is rotationally supported is robust to plausible variations 
in the gas geometry. 

4. COMPARISON WITH Ha KINEMATICS 

Both S0901 and the Clone have previously published kine- 
matic constraints in the form of Ha line widths (AuHa.obs)- To 
compare our [C ii] results to these, we use our kinematic models 
to predict the Ha line width Ai> h „ pre( i that would be expected, 
after accounting for the lower spectral resolutions of the pub- 
lished Ha observations. We do this by convolving our model 
line profile(s) with a Gaussian whose width <5 L L sf matches the 
reported line spread function (LSF) of the Ha observations, 
measuring the FWHM Aahglcohv of the convolved line pro- 
file, and calculating the LSF-corrected line width A7. H „ prcd = 
(Akjj a conv — <57L SF ) a5 . We finally convert the result back to veloc- 
ity units for ease of comparison (AuHospred = AkHa.pred x c/X n a ). 

For S0901, Hainline et al. (2009) report AuHa.obs = 308 ± 
12 kms -1 (FWHM), with a spectral resolution 5 A.lsf = 15 A 
FWHM, corresponding to 210kms -1 FWHM at 2.137 /xm 
wavelength. For our best-fitting model, the procedure described 
above yields Ai>Ha, P red = 266 km s -1 . If we apply the same 
procedure to the full set of models produced by the MCMC 
fitting, the mean and standard deviation of the result becomes 
Aupfa.pred = 255 ± 19 kms -1 . Combining in quadrature the 
scatter in our simulations and the uncertainty reported by 
Hainline et al. (2009), our models predict a line width that 
is smaller than the observed Ha line width at the 2a level. 

If the difference is real, there are a few possible explanations. 
First, the true LSF of the Keck spectrum could be 1 8-2 1 A 
(250-300 kms -1 ) FWHM, rather than the published measure- 
ment of 15 A (210 kms -1 ) FWHM (Hainline et al. 2009). 
Second, the [Cn] and Ha observations might sample differ- 
ent portions of the galactic disk. The Herschel observations 
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certainly sample the entire galaxy, while NIRSPEC’s (y.'76 slit 
could omit outer regions of the galactic disk. Third, the kinemat- 
ics of [C ii] and Ha could be genuinely different. If large velocity 
dispersions are associated with turbulence driven by star for- 
mation activity (Green et al. 2010; Lehnert et al. 2013), the 
kinematic effects of star formation would be most strongly 
apparent in Ha, which is closely associated with star forma- 
tion activity. The PDR gas that contributes strongly to [Cn] 
emission is found at the “skin” layers of molecular clouds 
and may yet retain the comparatively quiescent kinematics of 
molecular gas. 

Ha kinematic constraints for the Clone have been reported 
both by Hainline et al. (2009; slit spectroscopy) and Jones 
et al. (2010, 2013; integral field unit (IFU) spectroscopy). 
To compare to Hainline et al. (2009), we followed the same 
procedure we used for S0901. The best- fitting model gives an 
expected line width of Auna.pred = 160 km s' 1 FWHM, while 
the mean and standard deviation for all the MCMC runs was 
Auna.pred = 194 ± 52kms _1 . The observed Ha line width is 
consistent with these expectations, at AuHa.obs = 188 ±9 km s' 1 
FWHM (again with <U LSF = 15 A FWHM). 

The Ha IFU spectra reported by Jones et al. (2010, 2013) offer 
a more detailed comparison, with spatially resolved spectra at a 
spectral resolution of 52 km s' 1 . Their full measured rotation 
curve (i.e., the range of centroid velocities of Ha emission 
across the resolved image) spans Av = 159 ± 38kms _1 , with 
a line-of-sight velocity dispersion of 104 ± 37 km s 1 (Jones 
et al. 2013, Table 4). The results depend, again, on subtracting 
the instrumental LSF width in quadrature, but this time the 
LSF is narrower than the measured line width. The 104 km s' 1 
velocity dispersion corresponds directly to our a g parameter. 
The correspondence between the Jones et al. (2013) parameter 
Av = 159 ± 38 km s' 1 and our rotation speed parameter V \/2 is 
less direct. While we roughly expect Av Rs 2 i>i/ 2, our parameter 
refers to the rotation speed at the half-light radius of the [C ii] 
emission, while AV is measured across the full extent of the Ha 
emission. Still, we expect the kinematic parameters from Jones 
et al. (2013) to fall in the range 60kms _i < V \/2 < 80kms _1 
and o g Rs 104 km s' 1 . While this is far from our best- fitting 
model for the [Cii] rotation curve, it is within the 90% 
confidence region derived from the MCMC fitting results. The 
MCMC approach thus provides a better consistency check than 
a simple comparison of best-fit models. 

5. CONCLUSIONS 

Our [Cii] spectra of S0901 and the Clone showcase the 
potential of heterodyne spectroscopy for kinematic studies of 
high-redshift galaxies. Forster Schreiber et al. (2009) find that 
about 1 /3 of z ~ 2 galaxies are reasonably described as rotation- 
dominated disks, based on Ha IFU spectra. However, they are 
not able to resolve spectral features narrower than ~60 km s' 1 . 
Thus, our measured upper limits on the velocity dispersion 
of S0901, cs g < 23 km s' 1 (46 km s' 1 ) at la (2a), are one 
of the tightest limits to date on the velocity dispersion of a 
high-redshift galaxy. While our results are less definitive for 
the Clone, the maximum likelihood fits have a < 4 km s' 1 
(14km s' 1 ) at la (2a), and 65% of our simulations yield 
Av/ilog) Rs V\/ 2 /cr g > 0.4 — a (semi-empirical) criterion 
advocated by Forster Schreiber et al. (2009) for a disk to be 
considered rotationally supported. 

Spatially resolved observations at similar spectral resolu- 
tion and S/N will be possible for some <; > 2 galaxies us- 


ing (sub)millimeter interferometric line observations. This will 
enable cleaner measurements of the line-of-sight velocity dis- 
persion. ALMA will be the most powerful such instrument for 
the foreseeable future, offering spectral resolution comparable 
to HIFI, coupled with higher sensitivity and spatial resolution 
as fine as ~0’.'2. For our two particular objects, ALMA will not 
supersede our Herschel [Cii] observations: The [Cii] line in 
S0901 lies in a gap between ALMA’s frequency bands, and the 
Clone is too far north ( S = +51°) for ALMA observation. For- 
tunately, other strong lens systems will unlock the full potential 
of ALMA for this science. 

In conclusion, the Herschel observations of the [C ii] 158 /im 
line that we present here underscore the promise of the bright 
[Cii] line as a kinematic tracer for high-redshift galaxies. They 
provide a unique look at the internal dynamics of two z R* 2 
galaxies, and show that very small internal velocity dispersions 
can be found in the high-redshift universe. 
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